Reconstruction of the phase of matter-wave fields using a momentum resolved 

cross-correlation technique 
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We investigate the potential of the so-called XFROG cross-correlation technique originally devel- 
oped for ultrashort laser pulses for the recovery of the amplitude and phase of the condensate wave 
function of a Bose-Einstein condensate. Key features of the XFROG method are its high resolution, 
versatility and stability against noise and some sources of systematic errors. After showing how an 
analogue of XFROG can be realized for Bose-Einstein condensates, we illustrate its effectiveness 
in determining the amplitude and phase of the wave function of a vortex state. The impact of a 
reduction of the number of measurements and of typical sources of noise on the field reconstruction 
are also analyzed. 

PACS numbers: 39.20.+q,03.75.-b,06.90.+v 
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Many of the remarkable properties of atomic Bose- 
Einstein condensates (BECs) originate from the fact that 
in those systems a single wave-function is occupied by a 
macroscopic number of particles 0, 0, 0, EL 01 • The am- 
plitude of this complex wave function corresponds to the 
atomic density and can comparatively easily be measured 
using for example absorption or phase contrast imaging 
H E, II H 113 The phase of the wave function, on the 
other hand, is relatively hard to measure directly. Still, 
its measurement is essential for the characterization of 
many quantum-degenerate atomic systems, such as e.g. 
rotating BECs with vortices [Tlj . 

A situation in many ways similar is encountered in the 
context of ultrashort laser pulses of a few optical cycles 
duration. The characterization of these pulses, which 
are of considerable interest in both fundamental science 
and applications, requires likewise the knowledge of both 
amplitude and phase. 

The measurement of the time-dependent phase of ul- 
trashort laser pulses has found a solution that is in many 
respects optimal in the so-called Fre quen cy- Resolved Op- 
tical Gating (FROG) methods [H [H 111 E3 These 
methods, which can be adapted to many different situ- 
ations, offer very high resolution and precision, are sta- 
ble against noise, and can even compensate for or de- 
tect some sources of systematic errors. Single-shot mea- 
surements are straightforward fl6| , and measurements of 
fields with less than one photon per pulse on average 
have also been successfully demonstrated. Almost every 
aspect of the FROG-methods has been studied in depth. 
Good starting points for accessing the wealth of research 
literature on that subject are the review article 01 an d 
Ref. H3. 

The similarity of the phase-retrieval problems for ultra- 
short pulses and atomic condensates naturally leads one 
to ask whether the optical FROG techniques can be ex- 
tended to the atom-optical domain. This paper answers 
this question in the affirmative by demonstrating how 
a specific version of the general FROG method, called 
cross-correlation FROG (XFROG), can easily be trans- 
lated to the atom-optical case [ljj, |2(J, |2l| . 



Section H] briefly reviews the fundamentals of XFROG 
phase retrieval and discusses how it can be adapted to 
cold-atom scenarios. As an illustration, section ITU shows 
the reconstruction of the spatial phase of a rotating BEC 
with a vortex. We also investigate the sensitivity of the 
reconstruction to several sources of noise characteristic 
of cold atoms, as well as the impact of a reduction in the 
number of measurements on the field retrieval. 



I. PRINCIPLES OF THE XFROG METHOD 

The basic experimental setup for the field reconstruc- 
tion of an ultrashort laser pulse by means of the XFROG 
method is shown in Fig. ^ For simplicity we assume 
that only one direction of polarization of the electrical 
field needs to be considered. The unknown field is 
mixed with a known reference field if> T et(t) i n & nonlinear 
crystal with nonlinearity after a variable delay r. 
The fields i/j re f and ip should be of roughly comparable 
duration which means in practice that the pulse dura- 
tions can differ by up to about an order of magnitude. 
The sum frequency signal is 



(1) 



The form of that field shows that the reference pulse 
acts as a gate for the pulse, hence the acronym XFROG. 
The signal V'sig is then spectrally analyzed. The resulting 
spectrum, 



'XFROG 



dte- iwt ip(t)ip Iei (t 



(2) 



is the key quantity. It contains enough information to 
reconstruct amplitude and phase of the pulse tp(t). 

To see how this works |33| , we first note that it is suf- 
ficient to find ^/) s ig(t, r) since tp(t) can then be obtained 
by simply integrating over r — up to a multiplicative 
constant that can be determined from the normalization 
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FIG. 1: (Color online) Basic experimental setup for the 
XFROG method. G is a diffraction grating and the spectrum 
of the signal field is recorded with the CCD array. 



of ip. Thus the problem is cast in the form of a two- 
dimensional phase retrieval problem. Problems of this 
type have been studied for decades in the context of im- 
age restoration, see e.g. p3.l2^|. 

The field ip s i s can be found from Jxfrqg by means of 
the method of generalized projections [231. Y2M. l2fij. From 
Eqs. ||TJ and (J2J it is clear that V'sig must simultaneously 
belong to the two sets [34| 

O = {f(t, r)\f{t, t) = g(t)ip ie{ (t - t) for some g} (3) 

and 



F={f{t,r)\ 
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Hence, the solution is a field that satisfies the two con- 
straints 



^rig(t,r) eOnF. 



(5) 



Figure suggests that one can find the solution to 
this problem, called a feasibility problem in mathemat- 
ics and especially in optimization theory, by iteratively 
projecting onto the two constraint sets O and F. For 
closed convex sets with exactly one point of intersection, 
this method always leads to a unique solution. In our 
case, though, the constraint sets are not convex and their 
intersection consists of more than one function, as illus- 
trated in Fig. [2b- Hence, the algorithm of generalized 
projections is not guaranteed to converge for every ini- 
tial guess and even if it does the solution is not unique. 
Still, in practice it converges for a vast majority of initial 
guesses and the ambiguity in the solution is physically 
reasonable. (In particular, XFROG determines the field 
ipsig up to a constant phase.) On those rare occasions 
when the algorithm does not converge for a particular 
initial guess, this non-convergence is revealed by a large 
distance of the fixed point from the constraint sets in a 
sense that will be made precise below. Then one can sim- 
ply restart the reconstruction algorithm with a different 
initial guess. 

Specifically, ipsig is found from /xfrog a $ follows: 



FIG. 2: (Color online) Illustration of the constraint sets O 
and F (a) in the case where O and F are convex and their 
intersection contains exactly one element and (b) in the non- 
convex case with more than one element in the intersection. 
Fig. (a) also shows a typical trajectory of the generalized 
projections algorithm. 



1. Initialize (oj, r) with random numbers for its 
real and imaginary part. 

2. The projection onto the constraint set F is accom- 
plished by setting 
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which clearly guarantees that ip^ satisfies condi- 
tion (@J. 

3. Inverse Fourier transform tp^iuo^r) with respect 
to the first argument to find 



4. Determine ijj^ n \t) such that 
Z = I dtdr 



(7) 
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becomes a minimum. With ip^ n \t) determined this 
way form 



which is closer to the set O than ip^(t,r). 
5. Use the Fourier transform 



(9) 
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as a new input in step 2, and iterate until the error 
Z in Eq. JSJ) becomes sufficiently small. 

Upon exit from the algorithm, ij)( n > is the retrieved field. 
It turns out that in step 4 a simple line minimization 
along the gradient of Z with respect to is sufficient. 
A good starting point for that minimization is the mini- 
mizer of the previous step. 

An important characteristic of the XFROG method 
is its robustness, which results from two main reasons: 
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First, the XFROG signal contains a high degree of re- 
dundancy: If /xfrog is measured on a grid of size N xN 
the 27V unknowns of the field ip(t) - its real and imagi- 
nary parts - are retrieved from N 2 measured values. The 
XFROG algorithm makes use of this high degree of re- 
dundancy to yield a highly stable pulse retrieval. Second, 
the functional form of tp s i e is very restrictive in the sense 
that a randomly generated ^ s i g will not correspond to any 
physical pulse ip. Thus, if the XFROG signal has been 
deteriorated by systematic errors the XFROG algorithm 
will not converge for any initial guess and one can con- 
clude that the data is corrupted. Another advantage of 
XFROG is the extremely high temporal resolution that 
is achieved by making use of the Fourier domain infor- 
mation. Instead of being determined by the length of 
the reference pulse, it is essentially given by the response 
time of the non-linear medium. 

We now turn to the central point of this paper, which 
is to adapt the XFROG method to the characterization 
of matter-wave fields. We consider specifically the case 
of atomic bosons, with two internal states denoted by f 
and I. The atoms are assumed to be initially in a pure 
Bose-Einstein condensate (BEC) at temperature T = 0, 
with all atoms in internal state J.. The corresponding 
atomic field operator is 



V>(r) = cj.V(r), 



(11) 



where ip(r) is the condensate wave function and cj, is the 
bosonic annihilation operator for an atom in the conden- 
sate. 

The states f and j are coupled by a spatially dependent 
interaction of the generic form 

Vj = J d 3 rV(r — R)^>j(r)^(r) + H.c, (12) 

that is switched on at time t = 0. This could be provided 
for example by a two-photon Raman transition, with V 
then being proportional to the product of the mode func- 
tions of the two lasers driving the transition. 

For short enough times the ^-component of the atomic 
field is then 

f(r)txF(r-R)^(r). (13) 
The resulting momentum distribution of the f-atoms is 

n(q,R) oc y d 3 rd 3 r'e !;qr e~ 4qr ' ^j(r)^ T (r')} (14) 

x f d 3 rd 3 rV qr e- iqr ' 



Xi/j*(r)V*(r - R)V>(r')V(r' - R) (15) 



J d 3 re~ lclr i;(r)V(r - R) 



(16) 



When measured as a function of the shift R between cou- 
pling and condensate, the momentum distribution of the 
"["-atoms is therefore an XFROG signal, with the role of 



the reference field being played by the space dependent 
coupling strength. Hence, the XFROG algorithm can be 
used to fully recover the field ijj, n(q, R) being measured 
using for example absorption imaging after free expan- 
sion. 

In some sense, the situation for atoms is simpler than 
for ultrashort laser pulses. This is because a gate function 
of size comparable to the condensate size, such as e.g. the 
coupling strength V(r), is readily available. While in the 
case of photons the only possible gate is another ultra- 
short laser pulse, and hence one must use nonlinear mix- 
ing to obtain a signal suitable for XFROG, for atoms one 
can use processes that are linear in the atomic field. Fur- 
thermore, since the interactions between cold atoms tend 
to be stronger and offer a richer variety than is the case 
for light, a wide choice of processes can be employed to 
generate signals that can be used for FROG algorithms. 
Specific examples include s-wave interactions, the atom- 
optics analog of optical cubic nonlinearities, as well as 
the recently demonstrated coherent coupling of atoms to 
molecules, which corresponds to optical quadratic non- 
linearities m mm hi us H3 . 

A fully three-dimensional XFROG scheme as suggested 
by Eq. 116(1 is hard to realize in practice. First, it re- 
quires a very large number of measurements. Ten differ- 
ent shifts in each direction correspond to a total of 1000 
runs of the experiment. The large amount of data nec- 
essary for the fully three-dimensional scheme also poses 
serious challenges to the numerical reconstruction algo- 
rithm as far as computer memory and time are concerned. 
Second, in time-of-flight absorption imaging one typi- 
cally measures the column-integrated density so that the 
three-dimensional momentum distribution is not directly 
accessible. It appears therefore preferable to limit the 
reconstruction to have a two-dimensional scheme. 

It can be easily verified that if the field ip and the 
coupling strength can be factorized as 

^(r) = f(z)rl>(x, y), V(r) = h(z)V(x, y), (17) 

where z is the direction along which the imaging is done, 
Eq. (|16f) remains valid if we interpret r and R as two- 
dimensional vectors in the plane perpendicular to z and 
q as a corresponding two-dimensional momentum. Many 
fields of practical interest can at least approximately be 
written in the form Eq. i|17|) . F° r example, the two- 
photon Raman coupling mentioned earlier is of that type 
provided that the lasers are directed along z and the 
Rayleigh length is much longer than the extend of the 
atomic cloud in that direction. As a concrete example 
we demonstrate in the next section how the matter wave 
field of a rotating BEC with a single vortex can be 
directly reconstructed using the XFROG method. 



II. RECONSTRUCTION OF A VORTEX-FIELD 

We consider a BEC at temperature T = in a spheri- 
cal trap. We assume that the Thomas-Fermi approxima- 
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tion holds and that the healing length £ = (87rano)~ 1//2 , 
with a the s-wave scattering length and n the atomic 
density at the center of the trap, is much smaller than 
the Thomas-Fermi Radius R — (15A^ a a/a osc ) 1 / 5 a osc of 
the cloud. Here, N a is the number of atoms in the con- 
densate and dose is the oscillator length of the atoms in 
the trap. Under these assumptions the structure of the 
vortex core is essentially the same as that of a vortex 
in a uniform BEC and its wave function can to a good 
approximation be written as |3^ | 



rp(r,<p,z) = f{z) 



r/R 



(18) 

where we have used cylindrical coordinates with the vor- 
tex core at the symmetry axis. As discussed above, the 
^-dependence of the wave-function is unimportant and 
we will not regard it any further. In the following we 
use = 0.1, but none of our results depend strongly 
on this ratio as long as £/i£ <C 1. The real part of the 
wave-function (|18fl is shown in Fig. 

As an interaction Hamiltonian, or 'reference field' in 
the language of XFROG, we use Eq. (|12Jl with the Gaus- 
sian 



V(r) 



-{r/v 



(19) 



the ^-dependence being again irrelevant for our purposes. 

We have simulated an XFROG signal on a grid of 64 x 
64 points using ip of Eq. (|18fl and the reference field 
l|19jl. and applied the XFROG algorithm to reconstruct 
the field. In order the determine the optimal width of V 
we have repeated this procedure for various values of w. 
For each run we have determined the x 2 -error per degree 
of freedom |3fl 



1 



|V W (r«)-V(r«)| 



2N 2 ^ 

» >o 



,j grid points, 



(20) 

after 100 iterations. Here, N is the number of grid points 
in one direction. The results of these simulations are 
summarized in Fig. |21 They show that the XFROG al- 
gorithm works rather well for a wide range of widths pro- 
vided that they are comparable to the size of the conden- 
sate. As a rule of thumb, for our data \ 2 errors smaller 
than 10 -2 mean that the algorithm has qualitatively re- 
covered the original field. The comparatively poor qual- 
ity of the retrieved fields for larger widths is mainly due 
to unphysical correlations across the boundaries arising 
from the periodic boundary conditions that we are us- 
ing. These effects can in principle be avoided by using 
a larger grid. The best results were obtained for widths 
of w w 0.3R and in all that follows we use w = 0.35R. 
Larger widths tend to render the algorithm more stable, 
in the sense that it will converge to the correct solution 
for more initial guesses, while narrower reference fields 
result in faster convergence, but for fewer initial guesses. 

Several stages of the reconstruction algorithm are 
shown in figures EJb)-(d). (The field reconstruction in 




w/R 



FIG. 3: x 2 error of the reconstructed vortex field after 100 
iterations as a function of the width w of V. 



this example takes about 30 minutes on a Pentium 4 
CPU and uses approximately 400 MB of memory.) We 
show only the real part of the field, as the imaginary part 
shows a similar degree of agreement. The XFROG error 
Z of Eq. JHJ for the same simulation run is shown in 
Fig. [3] as a function of the number of iterations. Also 
shown is the deviation from the reconstructed field from 
the original field \ 2 ■ The figure shows that the algorithm 
converges exponentially after some initial stagnation. It 
also shows that, after the ambiguity in the total phase 
has been taken into account, the XFROG error Z is a 
good measure of the actual discrepancy between recon- 
structed and original field. This is important for real-life 
applications since in practice one does not know the orig- 
inal field so that x 2 cannot be calculated and only the 
XFROG error is accessible. 

Time of flight absorption imaging is a comparatively 
easy way to obtain many data points in the Fourier do- 
main, but it is cumbersome to obtain the data sets for 
different shifts R because each such set requires a new 
run of the experiment. Thus a natural question is the 
sensitivity of the algorithm to the number of displace- 
ments for which a measurement is performed. 

To answer this question we have calculated XFROG 
signals Ixfrog on grids of dimension n x n, interpolated 
them onto a larger grid of dimension N x N using cubic 
splines, and applied the XFROG algorithm to the result- 
ing XFROG signals. Fig. HJshows an example for n = 10. 
The recovered field shows good agreement with the orig- 
inal field, with differences in some details, e.g. near the 
maxima, resulting from the smoothing property of the 
interpolation with splines. 

To quantitatively characterize the dependence of the 
success of the field recovery on the number of measure- 
ments we have evaluated x 2 of the recovered field after 
100 iterations as a function of n. The result is shown 
in Fig. While the discrepancy between the recovered 
field and the original field grows as expected as the num- 




FIG. 4: (Color online) Reconstruction of a vortex state on a grid of 64 x 64 points. Figure (a) shows the real part of the original 
state and figures (b)-(d) show the initial guess and the reconstructed state after 50 and 100 iterations, respectively. 



ber of measurements decreases, the field is qualitatively 
retrieved down to n = 7. For even smaller n the field re- 
construction fails. In those cases the algorithm is honest 
enough to admit its failure by generating a large resid- 
ual XFROG error. Nonetheless, the conclusion is that 
a rather limited number of measurements is sufficient to 
accurately measure the field. 

We remark that these considerations are actually 
overly pessimistic. In reality one often knows the size of 
the atomic cloud and the XFROG signal is clearly zero if 
there is no overlap between coupling V and atomic field 
if). Thus, many data points in the XFROG signal can be 
padded with zeros. Even those zeros contain usable infor- 
mation for the XFROG algorithm because they encode a 
support constraint for the recovered field. Support con- 
straints, together with Fourier domain information, are 
widely used in image reconstruction for astronomic im- 
agery and under certain circumstances they characterize 
an image completely. Furthermore, among the remaining 
non-zero measurements some will give excessively small 



signals. In practice one can replace these measurements 
by other displacements that give a stronger signal. 

In an actual experiment the measurement of an 
XFROG signal will be subject to several sources of noise. 
In that case the intersection of the sets O and F in Fig. [21 
will in general be empty, and the XFROG algorithm de- 
termines a field that has the correct XFROG signal, i.e. 
belongs to the set F, and at the same time minimizes the 
distance to the set O as measured by Z. 

The XFROG signal Eq. ||SJ is insensitive to fluctua- 
tions in the total phase of the atomic field ip from one 
run of the experiment to the other, i.e. for different R 
— a consequence of the non-interferometric character of 
FROG methods in general. Hence it is sufficient to study 
the impact of uncertainties in the measurement of the dis- 
placements R themselves and in the total intensity of the 
XFROG signal from shot to shot. The latter can arise 
from variations of the number of atoms in the original 
condensate, from fluctuations in the coupling strength, 
and from fluctuations in the interaction times. We have 
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FIG. 5: FROG error Z and field error \ 2 f° r the simulation 
run of Fig. 0] 
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FIG. 8: (Color online) Same as Fig. E] but with simulated 
noise AR/R = 0.1 and AI/I = 0.1 in the XFROG signal. 
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FIG. 6: (Color online) Reconstruction of the field of Fig. [I] 
with the XFROG signal simulated on a grid of 10 x 10 and 
interpolated on a grid of 64 x 64 after 100 iterations. 




FIG. 7: x error of the recovered field after 100 iterations as 
a function of the number of measured data points n (same 
parameters as in figures 2] and |SJ . 



simulated XFROG signals with Gaussian fluctuations in 
R and in shot-to-shot total intensity / with variances AR 
and AI, respectively. Note that with our choice of nor- 
malization for ip and V the total intensity of the XFROG 
signals is of order one so that we can treat absolute errors 
and relative errors in I as being the same. 

We have applied the XFROG algorithm to these con- 
taminated signals and we have measured x 2 after 100 
iterations. To be more realistic, we have simulated the 
XFROG signals only on a grid of 10 x 10 and interpolated 
onto a grid of 64 x 64 as above. We have simulated the 
situation with errors in R only, errors in I only, and er- 
rors in both R and /. The results are summarized in Fig. 
El We were able to qualitatively reconstruct the atomic 
field up to relative errors as large as AR/R w 0.3 and 
AI/I = 1. The algorithm was found to be significantly 
less sensitive to uncertainties in the total intensity than to 
uncertainties in the displacements. Even XFROG signals 
with AI/I = 1 yield recovered fields of very high qual- 
ity. The relatively high residual errors for AI/I as 0.4 
and AI/I w 0.8 come about because the XFROG al- 
gorithm did not find the global minimum of Z for the 
particular initial guess. Restarting the algorithm with a 
different initial guess (but with the exact same contami- 
nated XFROG signal) lead to residual errors that nicely 
interpolate the the other data points. 

Figure |S| shows the retrieved field after 100 iterations 
for AR/R = 0.11 and AI/I = 0.11. The broad features 
of the original field are clearly reproduced and the noisy 
fine structure is exclusively due to the fluctuations in R 
and could be smoothed out by convolution with a Gaus- 
sian of width comparable to the healing length. Thus, 
although the error in the retrieved field increases with 
increasing noise, we conclude that the algorithm is not 
very sensitive to noise and yields at least reliable quali- 
tative information. In a sense, the field as recovered by 
the XFROG method contains much less noise than the 
input data. This is reminiscent of the situation encoun- 
tered in image reconstruction e.g. in astronomical appli- 
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FIG. 9: (Color online) \ error of the recovered field after 100 
iterations as a function of several noise sources. 

cations. And it was exactly for the purpose of removing 
noise from images by using a mixture of position space 
and Fourier space information that many of the methods 
described in this paper were first discussed. 

In summary we have shown how the powerful XFROG 



method of ultrafast optics can be adapted to phase mea- 
surement applications in ultra-cold atomic systems. In 
particular, this technique is capable of correctly recover- 
ing a vortex state and it is robust against noise, with a 
number of actual measurements that can be reduced to 
experimentally feasible values. 

An interesting question for the future is how the 
XFROG method can be used to study situations were 
there is no well defined phase of the atomic field for fun- 
damental reasons. For example the XFROG technique 
could be a valuable tool for diagnosing the transition of 
an elongated BEC (with long range phase coherence) to 
a quasi one-dimensional Tonks-Girardeau gas (without 
phase coherence), or the transition of a rapidly rotat- 
ing BEC into the quantum Hall regime. The XFROG 
method could also be of interest in the study of the so- 
called BCS-BEC crossover where a superfluid of Cooper 
pairs goes over into a BEC of tightly bound molecules. 

We thank H. Giessen for valuable criticism on the 
manuscript. This work has been supported in part by 
the US Office of Naval Research, by the National Science 
Foundation, by the US Army Research Office, and by the 
National Aeronautics and Space Administration. 
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